Spin-Orbit Coupling in an f-electron Tight-Binding Model 
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We extend a tight-binding method to include the effects of spin-orbit couphng, and apply it to 
the study of the electronic properties of the actinide elements Th, U, and Pu. These tight-binding 
parameters are determined for the fee crystal structure using the equivalent equilibrium volumes. 
In terms of the single particle energies and the electronic density of states, the overall quality of 
the tight-binding representation is excellent and of the same quality as without spin-orbit coupling. 
The values of the optimized tight-binding spin-orbit coupling parameters are comparable to those 
determined from purely atomic calculations. 

PACS numbers: 71.15.Ap, 71.15.Nc, 71.15.Rf, 71.20.Gj,71.70.Ej 
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I. INTRODUCTION 

The accurate determination of inter-atomic forces is 
crucial for almost all aspects of modeling the fundamen- 
tal behavior of materials. Whether one is interested in 
static equilibrium properties using Monte Carlo meth- 
ods, or time dependent phenomena using molecular dy- 
namics, the essential feature remains the origin, appli- 
cability, and transferability of the forces acting on the 
fundamental unit being modeled (atoms or molecules in 
most cases). First principles methods based on density 
functional theory have gained wide acceptance for their 
ease of use, relatively accurate determination of funda- 
mental properties, and high transferability. These tech- 
niques, however, are limited in their application by cur- 
rent computing technology to systems of a few hundred 
atoms or less (most commonly a few dozen atoms). Po- 
tentials that are classically derived (i.e., pair potentials) 
lack directional bonding (or at best add some bond angle 
information) and other quantum mechanical effects but 
are computationally far more tractable for larger simu- 
lations. Recent advances in tight-binding (TB) theory, 
which include directional bonding, but treat only the 
most important valence electrons shells, therefore show a 
great deal of promise. 

TB models have become a useful method for the 
computational modeling of materials properties thanks 
to their ability to incorporate quantum mechanics 
in a greatly simplified theoretical treatment, making 
large accurate simulations possible on modern digital 
computersiii^. Another advantage of these TB models 
is their ability to treat a general class of problems that 
include directional bonding between valence electrons, of 
particular importance for transition metal and /-electron 
materials. Finally, TB models are widely used in many- 
body formalisms for the one-electron part of the Hamilto- 
nian. It is therefore a useful representation of the band- 
structure for a more sophisticated treatment of electronic 



correlation, and has so been used-, for example, in dy- 
namical mean-field theory applications for Pu. 

In this report we present recent developments towards 
a transferable tight-binding total energy technique appli- 
cable to heavy metals. With the addition of spin-orbit 
coupling effects for angular momentum up to (and in- 
cluding) /-character, we demonstrate the applicability of 
this technique for the element Pu, of particular interest 
for its position near the half-filling point of the 5/ sub- 
shell in the actinide sequence and the boundary between 
localized and delocalized /-electrons^. 



II. TB METHOD 

The TB model used in this report is similar to that 
used in the handbook by Papaconstantopoulos'^. We 
have extended the calculations to include /-electrons^ 
and spin-orbit coupling^. As such, in this report we will 
elaborate only on those aspects of the technique that 
are unique to this work. A very brief recapitulation of 
the underlying TB method and its approximations is in- 
cluded to create the proper context for the addition of 
/-electrons and spin-orbit coupling. 

The Slater-Koster method^ consists of solving the sec- 
ular equation. 



(1) 



for the single-particle eigenvalues and orbitals, under the 
following restrictions: terms involving more than two 
centers are ignored, terms where the orbitals are on the 
same atomic site are taken as constants, and the result- 
ing reduced set of matrix elements are treated as variable 
parameters. The Hamiltonian, H , includes the labels for 
orbitals having generic quantum numbers a, (3 localized 
on atoms where the effective potential is assumed 
to be spherical, and can be represented as a sum over 
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atomic centers, 



site" terms, 



H. 
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where the on-site terms, Cq, represent terms in which two 
which we further decompose into "on-site" and "inter- orbitals share the same atomic site, and 

I 



E. 
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are the remaining energy integrals involving orbitals lo- 
cated on different atomic sites, and we have used transla- 
tional invariance to reduce the number of sums over bra- 
vais lattice points {Rn}, and the denote atomic basis 
vectors within the repeated lattice cells. Note that terms 
which have both orbitals located on the same site, but 
the effective potential (V^°*) on other sites have been ig- 
nored. These contributions are typically taken to be "en- 
vironmental" corrections to the on-site terms, and are not 
accounted for in the usual Slater-Koster formalism. For 
the inter-site terms, the two center approximation also 
consists of ignoring these additional terms in which the 
effective potential, V^^, does not lie on one of the atomic 
sites. Once this approximation has been made, the inter- 
atomic (i ^ j) matrix elements reduce to a simple sum 
over angular functions, Gwmi^ij), and functions which 
depend only upon the magnitude of the distances be- 
tween atoms, 



(5) 



where we have now adopted the usual convention of using 
the familiar Z,m angular momentum quantum numbers, 
and the axis connecting the atoms is the quantization 
axis. An equivalent expression for s^'m terms exists when 
non-orthogonal orbitals are used. The basis set used for 
the a and j3 quantum states are the cubic harmonics? 
whose functional forms are given in Table [J (with ap- 
propriate normalization factors) where |±) denotes the 
spin-state, which we will need for spin-orbit coupling. 

The Slater-Koster tables for the sp^d^ matrix elements 
can be found in standard referencesi^,, and we have used 
the tabulated results of Takegahara et al^ for the addi- 
tional matrix elements involving /-electrons. Typical TB 
applications are then reduced to using TB as an interpo- 
lation scheme; the matrix elements {hwm, swm and e^) 
are determined by fitting to ab-initio calculated quanti- 
ties such as the total energy and band energies. 

In this study we restrict ourselves to the determina- 
tion of optimal TB parameters at the neighbor distances 
in the face-centered cubic crystal structure (often used 
as a surrogate for the more complex ground state crystal 
structure of the actinides) near the equilibrium volume. 



Such tabulations have been extensively used^ in the study 
of materials with lower atomic number. To the best of 
our knowledge this is the first time that such parame- 
ters have been presented for light actinide elements that 
include the /-electron orbitals (although similar param- 
eters have been determined for the elements Ac and Th 
in an sp^d^ basis§). The TB parameter values so derived 
are available (on request) from the authors. 



A. Spin-orbit coupling 

The primary impact of spin-orbit coupling is to non- 
trivially couple electrons of different spin states, thus 
doubling the size of the TB Hamiltonian. The spin-orbit 
contribution to the Hamiltonian is given by 



(6) 



where ^(r) = {a'^ /{2r)){dV/dr), V is the total (crystal) 
potential. We neglect contributions from more than one 
center. A new Hamiltonian matrix can then be defined 
in terms of the spinless one. 



H = H + H' 



where 



H - \£,Lz 



(7) 



(8) 



is the spin-orbit coupling parameter between orbitals of 
orbital angular momentum I and primary quantum num- 
ber n located on the same atom, L± are the usual raising 
and lowering operators, and the azimuthal angular 
momentum operator, 



LzYlrr 



, (j)) = h^l{l + l)-m{m±l)Yi^±i 



hmYir 



The functions R'^i (r) are the non-relativistic radial wave 
functions. The spin-orbit contributions to the Hamilto- 
nian matrix can then be expressed in term of the TB 
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TABLE I: TB basis functions used for an sp d^f calculation. Note that /i(r) = 1/rK 



1=0 


1=1 




1=2 




1=3 


\s±) = v'l/47r|±) |pi±) = 


v/3/4,7/i(r)x|±) 


|di±> = 


\/5/l6TTf2(r)xy\±) 


l/i±) = 




b2±) = 


s/Vi^h(r)y\±) 


\d2±) = 


2V'l5/167r/2(r)a2:|±) 


\h±) = 


^7/WTTf3{r)x{bx^ - 3r^)|±) 


|P3±> = 




|d3±> = 


2v/l5/16Tr/2(r)za;|±) 


\h±) = 


v/7/167r/3(r)y(5y^ -3r^)|±) 






\di±) = 


v/l5/167r/2(r)(x"-a2)|±) 


\h±) = 


v/7/167r/3(r)2:(5z" -3r2)|±) 






M5±> = 


v/5/167r/2(r)(3z2 -r2)|±) 


1/5 ±) = 


v/l05/167r/3(r)a;(y2 _ 










l/6±) = 


v/l05/167r/3(r)y(z^ -2;^)|±) 










l/7±> = 


v/105/167r/3(r)^(2;=^ -a2)|±) 



basis functions listed in Table |T1 Rather than list contri- 
butions for the 32x32 matrix, here we list the matrices in 
the sub-blocks corresponding to each orbital angular mo- 
mentum. The p and d contributions have been previously 
discussed in relation to the tight-binding formalisir^iS^; 



to the best of our knowledge no / contribution has yet 
appeared in the literature. For completeness we detail 
the spin-orbit contribution for all values of the angular 
momentum up to / = 3. 
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where t = \/T5/2. B. Fitting tlie Parameters 

The values of the TB parameters were determined us- 
ing standard non-linear least squares optimization rou- 
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tines by matching energy band values derived from highly 
accurate first principles density functional theory (DFT) 
calculations^^. The technique is described in detail in 
a previous work^, where the DFT calculations in this 
case used a generalized gradient approximation DFT 
functionaU^, and the improved tetrahedron scheme^^ for 
Brillouin zone integrations. In this study we use as a 
starting point high quality fits to the scalar-relativistic 
energy bands and approximate atomic values of the spin- 
orbit parameters. The first step is to then use this fit for 
fitting the relativistic energy bands including spin-orbit 
coupling. Successive optimization steps then relax only 
the spin-orbit coupling paramaters (step 1), the remain- 
ing on-site parameters (step 2), and finally the inter-site 
terms (step 3). The fit quality through these steps is 
shown in Figure [TJ Note that the quality of the final fit 
is comparable to the original fit quality (open symbols 
at step 3) when only scalar-relativistic effects were taken 
into account. 




1 2 3 

Optimization Step 



FIG. 1: TB fit quality in terms of tlie cumulative root mean 
square (rms) errors at various steps of the optimization proce- 
dure. Step 1 relaxes the spin-orbit parametes {^„i), 2 relaxes 
the remaining on-site parameters, and 3 is a full relaxation of 
all parameters. Open symbols at Step 3 indicate the original 
scalar-relativistic fit quality. Note that the cumulative rms 
error is over all of the fitted bands (20 bands for Th, U, and 
Pu). Although the spin-orbit couphng is an atomic quantity, 
the improvement of our results in step 3 (which relaxes inter- 
site parameters) indicates some environmental effects should 
also be taken into account. 



calculations (the cumulative root mean square error in 
the TB fits to the first 20 energy bands in the irreducible 
Brillouin zone is 0.013, 0.013, and 0.072 Ry, respectively). 
Also note that we have included the "semi-core" 6p 
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III. APPLICATION TO THE LIGHT 
ACTINIDES, TH, U, AND PU 

A. Energy bands including spin-orbit coupling 

The first comparison between the TB fit and FLAPW 
calculations are the energy bands shown in Figure [2l 
Note the excellent agreement between the two sets of 



FIG. 2: TB energy bands for Th (a = 9.61), U (a = 8.22), and 
Pu (a = 8.14), shown in comparison with FLAPW valence 
energy bands (dotted lines). Note the excellent agreement. 
The abscissa for each calculation has been shifted such that 
the Fermi energy is at zero. Higher valence states (above the 
first 20) are not fit, hence the poorer fit quahty well above 
the Fermi level. 

states in the fit to better fix the available -p states in the 
TB basis. To expand the energy scale comparing the va- 
lence bands, the fit quality for the semi-core 6p states is 
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r A XZWQLATEKXSU 



FIG. 3: TB energy bands (dashed lines) for Pu semi-core 6p 
states, compared with FLAPW values (solid lines). 



shown separately in Figure [3] for Pu (all three elements 
have similar excellent fit quality for the more localized 6p 
states). Note that higher energy bands (well above the 
Fermi level) are not fit, hence the larger discrepancies for 
those levels. 



B. Density of states including spin-orbit coupling 

We also compare the total density of states (DOS) be- 
tween TB and FLAPW methods in Figure [H 

The TB method shown in the figure used a simple 
Fermi-Dirac temperature smearing method (with ksT ~ 
500) for integrating over the irreducible wedge of the 
Brillouin zone, while the FLAPW calculations used the 
improved tetrahedroni^ method with Gaussian smear- 
ing. From the comparison between the TB and FLAPW 
methods shown in the above figure, we note that the 
agreement is excellent, with all major features in the DOS 
reproduced by the TB calculations. There is a slight re- 
duction in the height of some of the larger peaks in the 
DOS for the TB technique, most likely due to the inabil- 
ity of the temperature smearing technique to represent 
the finer grained features as well as the improved tetra- 
hedron method does. 




e-E„ [eV] 



(a)Th 




(b)U 




C. Spin-orbit coupling terms 

It is interesting to compare the spin-orbit coupling pa- 
rameters, ^„i, predicted by TB theory for the various 
valence shells relative to the values predicted by accu- 
rate Hartree-Fock-Slater calculations of isolated atomsiS. 
This comparison is shown in Table [TTl 

Note the overall agreement between the TB fitted pa- 
rameters and the atomic values. The overall shift of a 
few tenths of an eV for the TB values is interesting, and 
this trend could be representative of crystal field effects 
(this speculation could be checked by performing equiv- 
alent fits at different densities). Equivalently, one can 



FIG. 4: TB (dotted lines) and FLAPW (solid lines) total 
DOS, including spin-orbit coupling. Note that the TB calcu- 
lation is in quite good agreement with the FLAPW results, 
despite using a different BZ integration method. The abscissa 
for each calculation has been shifted such that the Fermi en- 
ergy is at zero. 



compare the spin-orbit splitting of the electronic energy 
levels with the purely atomic case. This comparison is 
also shown in Table [TTl 
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TABLE II: Values of spin-orbit coupling strength, (^ni, and spin-orbit splittings, Ani = {21 + l)^„;/2, for the various valence 
electron shells predicted by the TB fit compared with purely atomic values using relativistic density functional theory (DFT)^, 
a Dirac-Slater atomic code (DIRAC)^*, and relativistic Hartree-Fock-Slater (HFS)^'' atomic calculations. Dashed entries are 
used for orbitals not populated in the atomic calculations. Values are in eV. 





Method 








Th 






As/ 






7 QA 


n on 




U.O-L 




U.DD 


DFT 


5.24 


7.86 


0.21 




0.52 






HFS 


4.09 


6.14 


0.30 




0.75 






TB 


4.19 


6.29 


0.20 


U 


0.51 


0.18 


0.62 


DIRAC 


5.96 


8.94 


0.19 




0.47 


0.24 


0.83 


DFT 


5.90 


8.85 


0.20 




0.50 


0.24 


0.84 


HFS 


4.38 


6.57 


0.30 




0.75 


0.35 


1.24 


TB 


4.64 


6.96 


0.23 


Pu 


0.58 


0.42 


1.48 


DIRAC 


6.92 


10.38 


0.20 




0.51 


0.31 


1.10 


DFT 
















HFS 


4.60 


6.90 








0.41 


1.43 


TB 


5.23 


7.84 


0.59 




1.46 


0.54 


1.90 





IV. CONCLUSIONS 

We have included /-electron and spin-orbit effects in 
a standard tight-binding method for solids in order to 
advance simpler simulation methods that are capable of 
the accuracy of more expensive, full-potential density- 
functional techniques. We have applied this TB tech- 
nique to elemental fee Th, U, and Pu, and have achieved 
excellent agreement with the electronic properties pre- 
dicted using a highly accurate FLAPW method. The 
fitted spin-orbit coupling parameters match very well 
the values independently predicted by atomic electronic 
structure calculations. This methodology bodes well for 
further TB investigations, especially for the study of de- 
fects, phonons, and dynamical properties. In future work 
we intend to develop a more transferable model based on 



a TB total energy formalism^, which should allow the 
straightforward calculation of detailed materials proper- 
ties. 
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